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Abstract 

We propose that reactor experiments could be used to constrain the environment dependence of neutrino mass 
and mixing parameters, which could be induced due to an acceleron coupling to matter fields. There are several 
short-baseline reactor experiment projects with different fractions of air and earth matter along the neutrino path. 
Moreover, the short baselines, in principle, allow the physical change of the material between source and detector. 
Hence, such experiments offer the possibility for a direct comparison of oscillations in air and matter. We demonstrate 
that for sin^(2Si3) > 0.04, two reactor experiments (one air, one matter) with baselines of at least 1.5 km can 
constrain any oscillation effect which is different in air and matter at the level of a few per cent. Furthermore, we 
find that using the same experiment while physically moving the material between source and detector improves 
systematics. 
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1. Introduction 

The concept of mass-varying neutrinos (MVNs) 
has been introduced by imposing a relation be- 
tween neutrinos and the dark energy of the Uni- 
verse [1-4] through a scalar field, the acceleron. 
Including the possibility of acceleron couplings 
to matter fields implies that the neutrino oscilla- 
tion parameters in vacuum and a medium could 
be very different [5], irrespective of the standard 
MSW efi'ect [6,7]. MVNs can have substantial im- 
plications for neutrino phenomenology, for exam- 
ple in the sun [8-10] or for cosmology and astro- 
physics [2, 11, 12], and they have been proposed as 
an explanation of the LSND anomaly [5, 13, 14]. 

Given existing data, precision tests of mass- 
varying effects are very difficult, since the different 
experiments are conducted under different condi- 
tions (different energies, baselines, matter distribu- 
tions, etc.). Practically, all evidence for neutrino 
oscillations so far involves neutrino paths in (so- 
lar or earth) matter, whereas no direct informa- 
tion is available on oscillation parameters in vac- 
uum or air. Moreover, the precise dependence of 
the neutrino masses on the matter density through 
the acceleron coupling is rather model dependent. 
The only plausible laboratory for a direct compar- 
ison between neutrino oscillations in different envi- 
ronments {e.g., matter and air) are short-baseline 



experiments. In this work, we consider the pos- 
sibility to constrain MVNs with reactor experi- 
ments, which are planned to measure the mixing 
angle On [15-17]. Such experiments are particu- 
larly suitable for this purpose, because one could 
think about building two very similar experiments 
with different material along the baseline, or even 
about moving the material along the baseline. In 
addition, reactor experiments suffer very little from 
correlations among the oscillation parameters, and 
new reactor experiments with identical near and 
far detectors will have an excellent sensitivity to 
sin^(20i3). Since the standard MSW effect is irrel- 
evant for such short baselines, any deviation of the 
energy spectrum between matter and air can then 
be interpreted in terms of MVNs. 

2. General formalism 

For the reactor experiments under consideration, 
the ordinary matter potential is, to a first ap- 
proximation, negligible because of the short base- 
line. The oscillation probability reads in the ex- 
pansion up to second order in sin(20i3) and a = 
AmlJAml, [cf, e.g., Ref. [18]) 



1 - Pee = sin2(26ii3) sin^ A31 + 
+ a''Al^ sin2(20i2). 



(1) 
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where A31 = Am'^iL/{4:E). Thus, one can easily 
see that short-baseline reactor experiments have 
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an excellent sensitivity to sin^(26'i3). In particu- 
lar, since o? ~ 10"'^, one can read ofF from this 
equation that for large sin^(20i3) ^ 10"'^ the sec- 
ond term is practically negligible, and we can use 
this formula in the two-flavor limit for analytical 
purposes. 

Suppose that we have any non-standard effect 
on neutrino oscillations on the Hamiltonian level 
which is different in air and matter. In this case, 
we have the usual vacuum Hamiltonian Tivac in 
air, and some non-standard Hamiltonian 7i,nat ~ 
'^vac + ^ns in matter. In general. Tins is a Her- 
mitian traceless'^ n x ti matrix (for details, see e.g. 
Ref. [19]). The Hamiltonian Timat can then be re- 
diagonalized in order to obtain the effective mixing 
angles ^13 and Am|]^ in matter. In order to pa- 
rameterize deviations between matter (tilde) and 
air (no tilde) parameters, we define 

8s = 013-^13- (2) 

For MVNs one has [8, 10] 

H„,at = ^ - ^i^P)\^ \^ - ^i^P)\ ' (3) 

where m = diag(mi) is the diagonal matrix of the 
neutrino masses and U is the mixing matrix, both 
in the background dominated environment (e.g., 
air), and M{p) is the mass matrix depending on the 
matter density p, which in general is non-diagonal 
in the basis used in Eq. Hence, one finds that 
in the case of MVNs Tivac and Tins have the same 
energy dependence, i.e., it is indeed the mass ma- 
trix which is modified. Note that according to 
Eq. Tins also depends on the m^. 

In the following we will assume that the neutri- 
nos propagate either through air or through mat- 
ter with constant density. This allows a model- 
independent test of MVNs, since only the density 
difference enters, and we need not to specify the 
detailed functional dependence of Mij{p). Fur- 
thermore, we assume that environment-dependent 
variations of Am|]^ and 612 can be neglected, and 
the two parameters and 5e suffice to describe 
MVN effects for reactor experiments under con- 
sideration. This assumption can be justified by re- 
quiring the consistency of solar neutrino and Kam- 
LAND data, due to the much higher matter densi- 
ties in the sun. The typical size of MVN parame- 
ters relevant for solar neutrinos [8, 10] is about one 
order of magnitude smaller than the sensitivity of 

''If it is not traceless, it can be made traceless by the sub- 
traction of a global phase. 



short-baseline reactor experiments, and in partic- 
ular, if the bounds derived in Ref. [10] are applied, 
MVN effects on Ani^i and 612 can be neglected 
(see also Sec. Ej). We stress that reactor experi- 
ments test different MVN parameters than solar 
neutrino experiments. 

Let us add, that for non-standard matter effects 
apart from MVNs (see Ref. [20] for a review) the 
energy dependence of Tins may be different than 
the one of Tivac- For instance, there may be a rel- 
ative factor 1/{2E) between Tivac and Tins leading 
to energy-dependent matter parameters similar to 
the standard MSW matter effect (c/., Ref. [19]). In 
this case, the assumption of an energy-independent 
effect can still be used as an approximation for 
short-baseline reactor experiments (by using the 
peak energy of the spectrum) as long as the neu- 
trino energy is far away from the MSW resonance 
energy. This requirement is only violated if the 
non-standard effects are more than two orders of 
magnitude larger than the standard MSW matter 
effect. 

3. Reactor experiment geometries 

In principle, we expect to obtain the best MVN 
sensitivities by the comparison of a matter to an 
air experiment. As we will see, using reactor exper- 
iments for the test of MVNs is mainly a matter of 
geometry. Some experiments have their baseline in 
air (such as Double Chooz; see, e.g., Refs. [13,14]), 
and others have it mostly in matter or mixed be- 
tween air and matter sections. In addition, mov- 
ing detectors have been discussed. Currently, three 
major principle schemes of reactor experiment lay- 
outs, which can be found in Fig.^ are: 

Type (a) Baseline mostly in air. Example: Dou- 
ble Chooz [21]. 

Type (b) Baseline entirely in matter because of 
flat terrain. Example: Braidwood [22], 
KASKA [23]. 

Type (c) Baseline mostly in matter because of 
hills. Examples: Daya Bay [24], Angra [25]. 

One could now think about three different concep- 
tual cases to test MVN effects. First, one could 
use existing (future) experiments with baselines in 
air and matter and combine them, such as Dou- 
ble Chooz or a larger reactor experiment in air 
combined with another reactor experiment in mat- 
ter, or the long-baseline beam experiments T2K 
or NOi^A. Second, one could modify the existing 
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(a) Mostly air 




(b) Flat terrain, underground 




(c) Hills, mostly underground 



Figure 1. Schematic illustration of different types 
of reactor experiment terrain layouts and potential 
modifications for the test of non-standard matter 
effects (not to scale). The filled boxes in (a) and (c) 
refer to additions of material in the far detector's 
line of sight, and the box in (b) refers to removing 
material in the far detector's line of sight after the 
primary operation period. 



experiments. In phase I, the experiment runs un- 
modified, and in a subsequent phase II the material 
in the line of sight of the far detector is changed. 
Examples of these modifications can be found in 
Fig.n For type (a), the far detector's line of sight 
could be covered with Earth or a rock wall. For 
type (b), a hole in the ground would lead to air 
propagation in the far detector's line of sight. And 
for type (c), the access tunnel to the far detector 
could a priori be built in the line of sight of far 
detector which is then covered with Earth or rock 
wall in phase II. Alternatively, one could drill an 
additional tunnel. As the last conceptual case, one 
could use a vertically or horizontally moving de- 
tector and the geometry of the site to change the 
material between source and detector. 

4. Analysis methods 

For the experiment simulation, we use a version 
of the GLoBES software [26] with some modifi- 
cations to simulate non-standard physics [27,28]. 
For the experiments, we mainly use the simula- 
tions of the Reactor-II and Double Chooz [21] se- 
tups from Refs. [16,29]. Reactor-II is an abstract 



background-free reactor experiment with identical 
near and far detectors, a baseline ifar — 1.7 km, 
and an integrated luminosity of 8000 1 GW yr. 
There are now many proposals for an experiment 
similar to this setup (c/., Ref. [17], and references 
therein). For Double Chooz, we use a background- 
free simulation with an integrated 60 000 unoscil- 
lated events, and a baseline Lfar = 1.05 km, which 
should be a very good approximation for the actual 
experiment. In addition, we will show some results 
in combination with NOi^A instead of a reactor ex- 
periment in matter, where we use the simulation 
based upon Refs. [30,31] adjusted to the proposal 
parameters [32] (L = 810 km, 12 km off-axis, 30 kt 
TASD) with a running time of five years in the 
neutrino mode. 

We simulate the data for a given set of "true" os- 
cillation parameters [33-36] Am§]^ = 2.2 •10~'^eV^, 
sin^ 26*23 = 1, Ato^i ^ 8.2 • 10~5eV^ sin^ 26*12 = 
0.83, and sin^(2^^l3) = 0.1 somewhat below the 
CHOOZ bound [37], Sqp — (relevant for beams 
only), and a normal hierarchy, and we assume no 
MVN effect in the "data", i.e., true 6a = — 0. 
Then we perform a fit to these simulated data, and 
we project the fit manifold onto the S^Se plane by 
marginalization of the standard oscillation param- 
eters. For this process, we impose external pre- 
cisions of 5% for each Am^i and 6*12 [34], which 
should be realistic estimates for the analysis time. 
We do not include the wrong hierarchy solution, 
because there is a good chance that the mass hi- 
erarchy will be measured on the timescale we are 
discussing (see, e.g., Refs. [29,32]).^ 

5. Systematics treatment 

Before we present our results, we need to discuss 
the reactor experiment systematics in greater de- 
tail. We use a treatment similar to Ref. [16] for 
a reactor experiment with two identical detectors, 
where the main systematics impact was identified 
as an effective normalization error of about 0.8%. 
Here we consider two different cases: 

1. One experiment with two phases: Phase I 
mainly in matter, phase II mainly in air. 
The detectors are not only identical, but also 
physically the same. 

2. Two different experiments, one mainly in air, 

■^Because reactor experiments are insensitive to the mass 
hierarchy, the results do, in principle, not depend on the 
mass hierarchy. However, the sign of the 5/\-efTects could 
change for a different mass hierarchy, because may come 
from a specific shift of one of the mass eigenstates. 
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the other mainly in matter. The detectors 
are identical, but not physically the same. 

We write the for both cases as 



E 



E 



[(l + 6 + 6,)iVf -O,] 



0^ 
b 

CTc 



(4) 



where Ni [Oi) are the oscillated (unoscillated) 
event numbers in energy bin i, the index x labels 
the data taken in air (a) and in matter (m), and 
one has to minimize Eq. with respect to 6, 6a, 
and bm- The error Uuc corresponds to an uncor- 
related normalization error between the air and 
matter data, whereas CTc is fully correlated. The 
interpretation of (T„c and depends on the setup 
and will be discussed below. 

Let us do some analytical estimations. First, 
we will see that the systematics mainly affects the 
s]x?{29\z) measurement, which means that we will 
be only interested in the s\r?{20i^) measurement 
in this section. Second, for simplicity, we con- 
sider only a rate measurement. Then the oscillated 
event rates can be written as N^' = (1 — S^fx)-^, 
where Sx = sm^{2df^), 



fx 



sm 



AmlL 



= 0{l), 



(5) 



and N = O \s the total number of events in one 
experiment. Using Sx,b,bx <C 1, the can be 
linearized. Since we are interested in the abso- 
lute error as on Sa — Sm, we then minimize the 
linearized Eq. (0J with respect to b, ba, and bm, 
and compute as taking into account the correla- 
tion matrix between Sa and Sm in order to obtain 



1 

77 



1 1 

J a J m 



1 1 

fa fm 



(6) 



For case 2 (two separate experiments), a^c cor- 
responds to the effective normalization error of 
each individual experiment, and the second term 
of Eq. © can be neglected because it is usually 
small compared to the first one. As one may ex- 
pect, in the (poor) statistics dominated regime 
^ cr„c, the error scales as l/V^, whereas 
for large event numbers l/V7\f ^ a^c the er- 
ror is limited by the systematical uncertainty auc 
{cf., also Refs. [38-40] for similar considerations). 



For Reactor-II, we have Af ^ 630 000 unoscil- 
lated events and a^c — 0.8%, which means that 
0.001 ^ auc, and we are in the system- 
atics dominated regime. 

For case 1 (one experiment), however, auc is ex- 
pected to be very small, since it only contains the 
time-dependent evolution of the systematics which 
is different for the near and far detectors. We 
only have to worry about values larger than about 
which is ~ 0.1% for a Rcactor-II setup, 
and we neglect it for the moment. In the statistics 
dominated regime l/\/N > (Tc, one can neglect 
the second term in Eq. ^ and one recovers the 
same limiting expression as for the two-experiment 
setup. For large M (and neglecting ct„c), the sec- 
ond term limits the accuracy. The scaling 
is cut off when the condition 



(/a frii) 

f 2 f2 
J a J 771 
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(7) 



is violated. Note that {fa — fm)"^ is a small 
number: Using Eq. one finds fa ~ fm — 
e{sm{AmlL/2E)), where e ^ (Ja/Atti^ < 0.1. 
Furthermore, also the second factor is small, since 
the fact that the experiments work at the oscilla- 
tion maximum implies A'm'^L/2E ~ tt. Using the 
numbers for Reactor-II, it is straight forward to 
estimate that the condition in Eq. is fulfilled, 
and hence, it is justified to neglect CTc for a one- 
experiment setup. We conclude that the case 1 
stays up to very high event numbers in the statis- 
tics dominated regime, and we simulate it in the 
next section by the most optimistic assumption of 
a systematics-free measurement. 

6. Phenomenological results 

Let us first discuss the case of using experiments 
proposed for measuring the standard oscillation 
parameters. We show in Fig. El left, the com- 
bined sensitivity of Double Chooz (air) and the re- 
actor experiment Reactor-II (matter). Obviously, 
there is a strong correlation between 6a and 60, 
which one can understand as follows: The exper- 
iment Reactor-II measures the matter parameters 
because of the much better statistics. For Sa > 0, 
the Am§i in air can be measured by Double Chooz 
by the spectral information. However, for Sa < 0, 
the oscillation phase becomes A31 <C 1, and the os- 
cillation probability in Eq. 0J can be expanded as 



sin^(20i3)Aii 



This implies that a 



smaller Sa (changing Am§]^) can be directly com- 
pensated by a larger Sg (changing ^13), which can 
be clearly seen in Fig. |21 left. The final precision 
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Reactor-ll+Double Chooz Reactor-ll+Reactor-ir (Twoexps) Reactor-ll+Reactor-ll* (Oneexp.) 
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i5i/(Am3, ),r i5A/(Am3, )„ i5A/(Am3, )„ 

Figure 2. The combined sensitivities to Sa and Sg (as relative changes between air and matter) for 
sin^(20i3)tr — 0.1 at the la, 2a, and 3a CL (2 d.o.f) for different experiments. The left panel refers to 
Double Chooz combined with a large reactor experiment Reactor-II in matter, the middle panel refers to 
two different experiments Reactor-II (matter) -|- Reactor-II* (100 m matter, 1500 m air, 100 m matter), 
and the right panel refers to the optimal case with the same experiment Reactor-II (matter, phase I) -|- 
Reactor-II* (phase II), i.e., actually the material is moved. 
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Figure 3. The sensitivities to 6a and 6g (as rel- 
ative changes between air and matter) at the 3a 
CL (1 d.o.f) as a function of the true sin^(20i3) 
for two different experiments Reactor-II (matter) 
-I- Reactor-II* (100 m matter, 1500 m air, 100 m 
matter) . 



for 6a and 60 is only comparable to the order of 
Am^i and sin^(20i3) themselves. 

As a next step, one could think about combin- 
ing two very similar reactor experiments, one in 
air and one almost in matter. We show in Fig. El 
middle, the results for this setup for two experi- 
ments with L = 1.7 km. Since both experiments 



practically have the same spectrum, the precision 
on 6a and 6g is extremely improved (below 10% 
at 3a, 1 d.o.f.). The correlation between and 
6e has almost completely disappeared, because of 
the long enough baseline. We find that there is 
no strong baseline dependence for 1.5 km < i < 
3.5 km, as can be seen from the lower part of Ta- 
ble n In Fig- ini we show the dependence of the 
sensitivity on sin^ (26*13). One observes that for 
sin^(26'i3) > 0.04, the difference of the parameters 
in air and matter can be constrained with a preci- 
sion of better than about 30% at the 3a CL.^ Sim- 
ilar results can be obtained for the combination of 
an air Double Chooz or Reactor-II experiment and 
a superbeam, such as T2K [41] or NOi^A [32] (c/.. 
Table In this case the sensitivity to 6g is some- 
what degraded because of additional correlations 
with 6cp and O23 in the beam experiment. 

Finally, we show in Fig.|2 right, the most opti- 
mistic results for the same experiment, i.e., using 
Reactor-II when 1.5 km of the matter is removed 
after the phase I run. For both phases we assume 
a running time of 5 years. As discussed in the pre- 
vious section, we simulate this situation by setting 
all systematical errors to zero, which leads to a 
further improvement of the sensitivity to 6g. As 
highlighted in Tabled such a setup allows a test 
of MVNs at the level of a few per cent at 3a CL. 



^Note that due to the small values of 613 ~ 0.1 a relative 
precision of 30% implies the impressive accuracy of ~ 0.03 
for Sg. 
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Experiment combination (5a/ {Arn^^)tr ^9/(^13)11 
Double Chooz (air) + Reactor-II (matter) 0.63 (0.29) 0.77 (0.29) 

Double Chooz (air) + ISlOvA (matter) 0.75 (0.30) > 1 (0.34) 

Reactor-II (matter) + Reactor-II* 0.09 (0.03) 0.12 (0.04) 
Reactor-II (matter), transformed into Reactor-II* 0.07 (0.02) 0.04 (0.01) 
Reactor-II (air) + 'NOvA (matter) 0.08 (0.03) 0.37 (0.23) 
Reactor-II -I- Reactor-II* , L = 1.0 km 0.33 (0.10) 0.20 (0.05) 
Reactor-II + Reactor-II*, L = 2.0km 0.08 (0.02) 0.13 (0.04) 
Reactor-II -I- Reactor-II*, L = 3.0km 0.07 (0.01) 0.13 (0.04) 

Table 1 

Relative sensitivities to 6a and 6g at the 3cr (Icr) CL (1 d.o.f.) for different experimental setups and 
sin^(20i3)tr = 0.1. If the upper and lower bounds are different, the maximum value is given. The definition 
of Reactor-II* is 100 m matter, 1500 m air, 100 m matter, whereas Reactor-II is in matter unless stated 
otherwise. The standard value for the baseline of Reactor-II is i = 1.7 km unless stated otherwise. 



7. Implications for mass-varying neu- 
trino models 

In the preceding section, we have demonstrated 
that reactor experiments could provide powerful 
tests of any neutrino oscillation model which is dif- 
ferent in air and matter in the Am3]^-sin^(2(?i3) sec- 
tor. Now we discuss some mass-varying neutrino 
model-dependent aspects. 

Let us first estimate the sensitivity to the el- 
ements M13, M33 of the environment-dependent 
mass matrix in Eq. For simplicity we take 

Mil = 0. Assuming an accuracy e ~ Ja/AttiIi ^ 
6e/0i3, a rough estimation gives a sensitivity 
of M13, M33 < e^/Arr^^. Now one can con- 
sider for example the following model for the 
environment-dependent mass matrix [8,10]: A/^ — 
X'^j X-^ n f / m'g , where / denotes a fermion in the 
medium, X'^j (X^) are the Yukawa couplings of 
neutrinos (the fermion /) to the acceleron, n/ is 
the number density of /, and mg is the acceleron 
mass. If we consider only the coupling to electrons, 
we have Ue — 6.4 • 10^ eV'^ in earth matter (for 
p = 2.8g/cm'^), and we find the following order 
of magnitude for the sensitivity of reactor experi- 
ments to the model parameters: 

This number has to be compared to the value 
3 • 10~^® obtained in Ref. [10] for the sensitivity of 
solar neutrino data to MVN parameters in the 1-2 
sector. Naively one expects that solar neutrino ex- 
periments are much more sensitive to MVN effects 
than reactor experiments, since there is roughly 
a factor of 100 between the densities in the solar 
core and earth matter, and there is another fac- 
tor -y/Amgj^/Aml]^ ~ 5, which characterizes the 



mass scale sensitivity of the experiments. However, 
these factors are partially compensated by the high 
precision of reactor experiments, e ~ 0.03. 

In principle solar -l-KamLAND neutrino data 
have some sensitivity to ^13 [33-36]. The above 
argument on the factor 100 between the densi- 
ties in the sun and the earth suggests that solar 
data should also be able to provide some informa- 
tion on MVN effects for ^13. Such a constraint 
requires model-dependent assumptions about the 
density dependence, since one has to relate the 
environment-dependent values of ^13 inside the 
sun, in earth matter, and in air. Depending 
on such model assumptions, one could have con- 
straints on 6g from solar neutrino data, which could 
change the interpretation of our results: For 6e 
fixed to 0, already planned experiments in mat- 
ter in combination with Double Chooz will provide 
very good results for 6 a {cf-, left panel of Fig. |2J). 

A very particular mass-varying neutrino model 
has been considered in Ref. [14] to explain the 
LSND experiment and a potential null result in 
MiniBOONE by the combination of mass-varying 
neutrinos and one additional sterile neutrino mix- 
ing with the active ones in matter. Since this 
model predicts no oscillations in air for short- 
baseline reactor experiments, a signal for sin^ (2^13) 
in air (for example in Double Chooz) would be 
a strong rejection of this model, i.e., 6a has to 
be tested only of the order of one. Furthermore, 
our setups can test any additional sterile neutrinos 
with mixings different in air and matter [11,14] as 
follows: If the fast oscillation frequency associated 
with the sterile neutrinos AM^ leads to oscillation 
effects already at the near detector site (such as 
for AM^ > 1 eV^), then both the near and far de- 
tectors will observe a reduced overall flux. This 
can neither be faked by sin^(26'i3) nor Amfi being 
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different in matter. Thus, a sterile neutrino con- 
tribution in matter can be measured on the level 
the reactor flux is known, typically some per cent, 
see e.g., Ref. [42]. 

8. Summary and conclusions 

We have considered several reactor experiment se- 
tups with baselines in air and matter to constrain 
any non-standard contribution to the Hamiltonian 
which is different in air and matter. In particular, 
we find that new reactor experiments with near 
and far detectors could provide stringent bounds 
for mass- varying neutrino models, which lead to 
environment dependent effects in the 1-3 sector 
(Ato§i or sin^(26'i3)). 

Non-trivial constraints can be obtained already 
from Double Chooz (mainly air) combined with a 
matter experiment, although in this case correla- 
tions between the parameters limit the sensitivity. 
Best results, however, are obtained if two identical 
setups with baselines of at least 1.5 km with sub- 
stantially different material between the near and 
far detectors are compared. For sin^(20i3) > 0.04, 
deviations of oscillation parameters in matter and 
air can be constrained at the level of few per cent. 
Using the same experiment and physically moving 
the matter between the two detectors after the ini- 
tial operation period decreases the impact of sys- 
tematical errors. In this case, the relevant sys- 
tematics issue is the uncorrelated time-dependent 
change between the near and far detectors, which 
should be extremely small for practical purposes. 
We conclude that new reactor experiments could 
be excellent candidates for the test of mass- varying 
neutrinos provided that sin^(26'i3) is not to small. 
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